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ABSTRACT 

Modified Laplace transformation method is applied to component (f)^ theory and the 
finite temperature problem in the massless limit is re-examined in the large N limit. We 
perform perturbation expansion of the dressed thermal mass in the massive case to several 
orders and try the massless approximation with the help of modified Laplace transforma- 
tion. The contribution with fractional power of the coupling constant is recovered from 
the truncated massive series. The use of inverse Laplace transformation with respect to 
the mass square is crucial in evaluating the coefficients of fractional power terms. 
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1 Introduction 



It is well known that the conventional perturbation theory for bosons breaks down when 
the temperature is sufficiently high compared to the mass scale of the theory under 
consideration^ . The origin is easily seen in the massless limit as the emergence of infra-red 
divergences of Feynman diagrams. This is linked to the fact that the perturbative infra-red 
structure becomes worse at finite temperature. For example, the boson propagator shows 
the behavior more singular at small momentum region than at zero temperature, and the 
temperature correction remains finite for the self-energy while at zero temperature it van- 
ishes at zero mass in dimensional regularization. The nonvanishing of mass correction 
activates the hidden infra-red divergence of sub-diagrams coupled to the tadpole, and thus 
the divergence manifests itself in Feynman amplitudes. 

The presence of the thermal mass means that there would be no infra-red problem if one 
can incorporate such effect from the outset of perturbation expansion. It also means, as 
well known, that the ordinary perturbation expansion in powers of the coupling constant 
breaks down. For example, in a single component 0^ theory with interaction, the 
pressure P is given at temperature T (= by 

-TT^ 1 1 



90 48 127r ^ ' 



(1) 



The presence of the fractional power of the coupling constant A exhibits the break down 
of perturbation expansion. 

To obtain the amplitude free from the infra-red singularity, one may sum up a set of 
Feynman diagrams to all orders. Similar but more systematic approach would be to shift 
the perturbative vacuum such that the fields are thermally screened at the tree leveP'^'^'^. 
This program is implemented by adding the thermal mass to be induced by the loop 
correction to the free part and subtracting from the interaction part. Although this is 
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physically appealing, it is not effective in the study of magnetic part of gluon self-energy, 
because the one-loop contribution is absent and higher order diagrams are plagued with 
the infra-red singularity. This is a reason why we look for another expansion scheme. 

Recently a new perturbative scheme was proposed^'^ and it seems to be a convenient 
scheme to deal with the infra-red problem at finite temperature. The scheme is based on the 
use of the Laplace transformation which is modified as to fit the perturbative framework. 
In the case of anharmonic oscillator, the infra-red problem at zero temperature is well 
controlled by the modified Laplace transformation (MLT) method^. The purpose of the 
present paper is to apply the method to the finite temperature perturbation theory of N 
component theory at D = A ^. To check our result we confine ourselves with the 
leading order oil/N expansion. We consider the self-energy and investigate how the non- 
analytic structure in the coupling constant similar to (1) can be recovered from massive 
perturbation expansion which gives a simple power series in the coupling constant. We 
do not sum up Feynman diagrams to all orders, but just deal with the truncated series. 
Nevertheless we will show that first several coefficients of the fractional power terms can 
be obtained from the truncated perturbative series. 

2 Survey of the modified Laplace transformation method 

For the self-containedness we here survey the basic points of MLT method^'^. Since we 
confine ourselves with the approximation of a massless theory, our presentation emphasizes 
the relevant respects only. 

Let /((t) be a physical function and a be monotonic with respect to the mass. Then 
the massless limit is equivalent with the limit, cr — > 0. Suppose that /(a) has perturbative 
expansion in l/cr and is given at the perturbative order k a^s fk{a) — YLn=o ^n/c*^ where 
p is some constant. Then the "infra-red problem" occurs when we ask the behavior of f{a) 
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for small a. This is the situation we encounter at finite temperature theory. Since we 
cannot take the massless limit in /^(cr), we may fix cr at some small value where the lower 
limit of perturbation series is signaled. In this approximation scheme, it is convenient to 
consider the inverse Laplace or Heaviside transform with respect to d as we can see below. 
Heaviside transform of /(cr) is given by the Bromwich integral, 

/•s+ioo ^(j \ 

M= —exp{ax)-f{a), (2) 

Js-ioo Zm (7 

where the parameter s represents the location of the vertical contour. The contour of 
integration should be placed on the right side to all the possible poles and the cut of 
f{a)/a. Then, if x < 0, the contour may be closed into the right half circle and f{x) is 
found to vanish. As well known, f{x) is the kernel of Laplace transformation of /(cr). We 
have 

/oo 
dxexp{—ax)f{x). (3) 
-oo 

Since f{x) — when x < 0, the integration range reduces to [0, oo). However f{x) involves 
the step and S functions and it is generally convenient to keep the range as (— oo, +oo) to 
handle integration easily. 

In our massless approximation, we begin with the observation that 

Ji^o/(^) = Jl^/>)> (4) 

where we assumed the existence of the limits. This ensures us that /(x) can be equally used 
for the massless approximation. Let us denote the Heaviside transform of /^(cr) as fk{x). 
Although we cannot take the massless limit for /^(x), we may approximate /(oo) by fixing 
X some large value, x*, in fk{x). Then we note that, in some cases including our finite 
temperature problem, the transformed function has an advantage that the convergence 
radius is larger than /(cr). For example, when / = Z^On/cr", / = X^(a„/n!)x". Therefore, 



4 



if lima;^oo f{x) converges, then the transformed function gives us more chance to know the 
massless value by increasing the order of expansion. 

To summarize, calculating fk{x) from (2), we approximate /(O) by fk{x*). The exphcit 
way of fixing x* will be discussed in the next section. 

This approach can be naturally extended to the calculation of the massive case''. Let 
fk{a,x*) be the A;-th order approximation of /(cr). Then fk{a,x*) is given by 

fk{cr,x*) ^e:icp{-ax*)fk{x*) + a dxex.p{-ax)fk{x). (5) 

J —CO 

The right-hand-side defines our modification of Laplace transformation. In addition to the 
naive cut-off term (the second term), we are led to add the first term. This term is critical 
in approximating /(a) for small a. 

In performing Heaviside transformation, the following formula is useful; 

where 6{x) denotes the step function {6{x) = for x < and 1 for x > 0). For larger p 
one can obtain the transform by using 

For example we have 

^ (n = l,2,3,---) 



In this way we find 



where we have omitted the 9 and 5 functions. We note that for fractional p the transform 
does not vanish. 
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3 Application of MLT method 

We work with the imaginary time formahsm of finite temperature field theory. The model 
we consider is defined by®: 

= \{dct>f + ^mV^ + + Ce,c, (10) 

where 0^ = Y.n=i 0n ^"^^ ^e,c denotes the collection of counter terms. We use dimensional 
regularization at L> = 4 — 2e with MS prescription. Before turning to detailed calculation, 
we spend some words on our basic strategy. 

We write the massive perturbation expansion of the thermal mass to order k as 

k 

^rri^ + Y, Cn{m)g'^. (11) 

n=l 

Since we are interested in the small mass behavior of massive perturbation theory, we 
expand each loop correction in the mass. This is equivalent to use the high temperature 
expansion. For example, the one-loop contribution to the thermal mass gives 
8C(2) 47rm , Stt^/^ ^ (-1)" ./3Wn 



Cl = g 

= 9 
where 



8C(2) 47rm , ^ ,3, 



/32 /3 



12) 



The term linear in m is the crucial one in this correction. It is in a sense non-perturbative, 
since the mass enters in the Lagrangian as the square of the mass, m^. Actually all the 
infra-red divergent behavior of sub-diagrams with many legs comes from that term by 
the differentiation with respect to rn?. 

There is a subtlety in choosing the transformation variable cr. For instance we can 
enlarge the convergence radius of the thermal mass, M^(m), in both cases where a — m 



and u = m^. In this paper, we choose the choice, a — m^, according to the successful 
choice in the anharmonic osciUator^. 

Heaviside transform of Ci (m) , Ci (x) , is then given by 

-8C(2) 47r 



ci{x) = g 



(14) 



/5v/ir(l/2) 

where Reven{x) denotes the Heaviside transform of Reveni'm^)- Since any term in Reveni'm'^) 
is an even positive power of m, any finite sum involved gives zero (see (9)). Although this 
does not necessarily mean that Reven{x) is zero, we may safely neglect it when we focus on 
the massless approximation. On the other hand, while the second term, 47r//9Y^r(l/2), 
tends to zero when x — > oo, we keep it because it gives non-negligible contribution at finite 
perturbative order where x is fixed to some finite value in the massless approximation. To 
summarize, we keep terms as the transformation suggests: Terms needed are of negative, 
zero and odd powers of m (which corresponds to positive, zero and fractional powers of x 

)• 

Before performing our approximate calculation to the full several loops, let us illustrate 
in the next subsection how the leading fractional term of order g^^^ comes out in our 
approach. 

3.1 Calculation of g^^"^ term 

It is well known that the leading fractional term of order g^^^ comes from the series of 
Feynman diagrams shown in Fig.l, each of which gives most dominant contribution for 
small m at each order in g. Let us first focus on the dominant graphs and try to approximate 
the leading fractional term. 

The diagrams shown in Fig.l gives: 



8C(2) 



(5^ 



3(3m^ 27r , 8C(2)7r g^ (8C(2))27r 1 ■ 3^, 



i{l-^)--—g^ + 



TT ' iijxi (3^m^ 2! /3^m^ 3!2 
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(8C(2))37rl-3-5 , 



+ (sub dominant terms), 



(15) 



(3^m^ 4!22 

where (sub dominant terms) denotes the sub-dominant pieces in the small m expansion. 
Hereafter we omit those sub-dominant terms. We have kept the next-to-the leading term 
(the second one in (15)) in the one-loop contribution. This is because the leading con- 
tribution of each higher order diagram comes from this term and the power counting of 
higher terms with respect to m suggests that it is the first term in the dominant series. 
From (6) the Heaviside transform of M^^^-^^^^ with respect to is found to give 



■^^■^ dominant ) 



K(2) 



3/i 



V^r(i/2)^ /?r(3/2)^ ' r(5/2)/33 2! 



9' + 



/32 

_ (8C(2))W/M -3 
r(7/2)/35 3!2^ 

Nontrivial result is obtained from the 2-loop order. At 2-loop order, we have 



(16) 



■^dominanti-^) 



802) 
/52 



.(1- 



3/3 



2ti\/x 9 



(17) 



7rV^r(l/2)^ r(3/2)/5^ 
The function M^^^^^^^^{x) behaves as shown in Fig.2 and the limitation of the result is 
seen, for example for g — 0.05, at ~ 1.52. When x/P'^ exceeds 1.52, the highest term 
in (17) dominates over and therefore the result (17) loses its validity for larger x. Thus, 
the upper limit of perturbative region may be given at the stationary point of -/Wj^^j^^^^ {x) 
and it is specified by 



^Afi_,(a:) 8C(2) 



dx 



3gP 



Tvg 



0. 



(18) 



.27rr(l/2)x3/2 r(3/2)/?,TV2 
We note that the first term in (17), which is independent of x, does not enter the station 
arity condition (18). The solution of (18) is given by 

3/32 



X 



4:Tr^g' 



(19) 



Note that x* is proportional to 51 ^. Hence the fractional power in x produces the fractional 
power of g. This is the mechanism to generate the fractional power of g which generally 



holds in our approach. By substituting x* into M^^^^^^^^^x) we obtain the approximation 
of Mdominanti'n^ = 0)- The result reads 



-^^-^ dominant K"^ I 



8C(2) 

/32 



9 - 



4v^ 
0r 



3/2 



(20) 



Thus, we have recovered the g^^'^ term from the truncated power series in g. Since the 
exact result reads M'^{m = 0) = %i{2)li?\g - 2y/3g'-^/^ + 0{g^)], the first term in (20) is 
exact and the error in 0{g^/'^) is about 13 percents. This is satisfactory as the lowest order 
approximation. 

Now we investigate the effect of multi-loop dominant graphs. Since the stationarity 
condition always gives x* as (3"^ / {C,{2)g) times the number, it is convenient to define 

and consider the transformation with respect to a. We note that m^/g is renormalization 
group (RG) invariant and therefore a is. In terms of a, (15) is written as 

8C(2) 



dominant 



g - VSg'^'Aia) 



22) 



where 



1 11 11-3 11-3-5 
^ ^ ~ ^2!2 ^ ^3!22 ~ ^ 4!23 ^ 



■ - (-1)" (2n-l)!! 

^^2^ „n- 



n=0 



cT"+i (n+ 1)!2" 



Let us truncate (23) to the k + 2 loop level. The truncated series is given by 

A(n\ 0./3I f (-ir (2n-l)!! 

The corresponding Heaviside function Ak{x) is given by 

^'^^^^ r(l/2)xV2 ^ ^0 (n + l)!2-r(n + 3/2)"^ " 



(23) 



(24) 



(25) 



Here x is the conjugate variable of cr defined by (21) and should not be taken as the previous 
X appeared in (16), (17), (18) and (19). We emphasize that, while the convergence radius 



of A{a) is 1, that of A{x) is infinite. As mentioned in the introduction, this is the essential 
advantage in deahng with Ak{x) rather than Ak{a). Our task is to approximate ^(0), which 
is equal with 2, by Ak{x) where x may be optimized by the stationarity requirement. The 
stationarity condition to fix x leads to dAk{x)/dx = and the solution is found at even 
orders. By substituting the solution into Ak{x), we obtain the following approximation for 
A{0) = 2; 



2 


— loop 


2.25676 


[x 


= 1), 


4 


— loop 


2.06054 


[x 


= 1.596), 


6 


— loop 


2.02096 


[x 


= 2.181), 


8 


— loop 


2.00835 


[x 


= 2.759), 


10 


— loop 


2.00358 


[x 


= 3.334). 



(26) 

Thus, the sequence gives the good approximation of ^(0) and gf^/^ term. Finally we point 
out that the omitting of the first term in A}^{a) (see (24)) gives poor result and its necessity 
is confirmed. 

3.2 Approximation via full loop expansion to 10-loops 

The previous procedure can be applied to the full loop corrections of the dressed thermal 
mass. 

We first sum up the logarithmic corrections by fixing the value of t. Let us fix t simply 

by 

i = 7 + log^ = 0. (27) 
This gives //^ = 47re~'''/ and the coupling constant g is given by 



9 = 



log ^ + 7^ 



(28) 
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where A denotes the intrinsic scale in MS scheme. If one wishes to recover the ordinary 
logarithmic terms, one can do it by expanding gfj as gp = g{jji) {l + tg{jji) + t^g^{jji) + •••). 
In this case, RG invariance will be satisfied to the same level of ordinary perturbative 
expansion. 

From 2-loop level, the terms with odd power of m appear to all orders in the high 
temperature expansion. For instance 



(29) 



8C(2)- " 3^ " ^6^^^^'^ + 128^^^'^ + • • • J 

where ■ ■ ■ denotes the higher order terms in m and even powers of m is fully omitted. The 
terms with positive and odd powers of m survive under the Heaviside transformation with 
respect to m^. Hence, although those terms go to zero as cr — > 0, they are needed at finite 
order masslcss approximation. Wc note that any term in the small mass expansion of C2 
is of ,|-order, but the higher order term in m in fact contributes to the higher order in gp 
of massless approximant as suggested in (29). 

In small mass or high-temperature expansion, m independent pieces (the first term in 
(12) and (29), for example) play a special role under our approach because those are not 
affected by our approximation scheme and go through to the final result. Those terms are 
collected as 

8C(2) 



5/3 + 65| + ^C(%^ + 0(4) 



= (30) 



(3^ 

and give the purely power like part of M^(m = 0). Actually the above result agrees with 
the exact result. On the other hand, the terms with fractional power of cr give rise to 
the fractional power (in g) part of the approximant of M^{m — 0). At this stage, let us 
summarize the classification of terms responsible for the approximation of M'^ijn — 0). 
The terms zero-th order in small mass expansion were gathered into M| and exactly gives 

11 



the integer power part of M^{m = 0). The terms of odd power of m are gathered into 
and contribute to the fractional power part of M'^{m = 0). Thus, neglecting terms of 
even positive powers of m (including the tree-level contribution), we can put 

M2(m) = Ml + M^{m). (31) 

Now we explicitly carry out the approximation of the massless value of M"^. We have 
calculated the radiative corrections c„ up to 10- loops and classified the result in powers of 
gi3. Then we can write 



A{a) + gf,B{a) + gp{a) + 0{g^ 



(32) 



We further classified coefficient functions, A{a), B{a), ■ ■ ■, based on ( functions. The result 
is shown in appendix. We thus arrive at 



M'^ = -V34'=[2« + 369,+ (5|ac(3)-5;)j|+(^C(3)-I|5>C(5) + |<i)jJ 
+ (^C(3) + %lic(3)^ - ^C(5) + %iC(7) - + 0(4)]m 



6 ' 4/^^ V 4 ^ 36 ' 8 

,e - ,-(3)-^C(5) + ^C(r)-l^ 

Our massless approximation goes as in the previous subsection: We first perform the 
Heaviside transformation of M^(m) which gives M^{x) = + M'^{x). Here M'^{x) is 
given by the series same as (33) except that the coefficient functions a{a), b{a), ■ ■ ■ are 
now replaced by the Heaviside transforms, a{x), b{x), ■ ■ ■. Then we look for the stationary 
points of them to evaluate their massless values. The result is shown in Table 1. 

We find that, at the 10-loop level, the first 5 coefficient functions. A, B, C, D and E are 
approximately evaluated. The coefficients a, 6, c^, and already reach to the sufficient 
accuracy at this level. The accuracy of c, d and e is however not so good. The reason is 
that c{x),d{x) and e{x) start with x^l'^^x^^'^ and x^l'^ terms, respectively, and therefore 
need higher order contributions to roughly reach their values at x = oo. Prom Table 1, we 
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obtain the following appi oximant of the fractional power part of (m = 0) , 



2.0036 + 2.9220t/^ - 0.6172c/| + 11.2082^| + 8.7295^1 + 0(c/| 



(34) 



The exact result reads, 



Mj^(m = 0) = -V34/' 



2 + ?,gi3 - 0.8075c/| + 13.6914^ + 4.0193^^ + 0(c/|) 



(35) 



In view of (34) and (35), we can say that our approach can produce the fractional power 
part of the self-energy to the satisfactory level. 

4 Conclusion 

In the large N limit of N component 0^ theory, we found the mechanism that how the 
integer and the fractional power terms in the massless theory are systematically classified 
and generated via massive perturbation theory. The fractional terms comes from the odd 
power terms of m in the high temperature expansion. To 10-loop order, the terms of order 
g^J"^ , g^J"^ ^ gy^ , g^ and g^^"^ were approximately calculated with the help of MLT method. 
Although only the two point function is considered in this paper, our approach would be 
directly applied to the thermodynamic functions, without referring to the mass correction. 

In our approach, we can deal with all the fractional power terms on equal footing. 
We have not performed infinite sum which will be difficult in general. The characteristic 
features in our approach are that the evaluation of the massless values is based on the 
truncated sum which can be obtained for general theories and, even when the one-loop 
contribution vanishes at m = 0, our evaluation procedure works. Therefore we expect 
that our approach may be effective for more complicated theories. Long time ago, Linde 
observed that perturbative hot QCD breaks down at and above a certain perturbative 
order^, since the magnetic gluon mass would be of order g^T and then the infinite number 
of diagrams contribute to a single order in g (See also refs. 10,1). Our approach may 
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provide a new calculational framework toward the resolution of the problem. Although 
the calculation of higher order diagrams is tedious and involved, the fundamental problems 
issued by Linde become harmless. 



Appendix 

In this appendix we summarize our perturbative calculation. The function / is represented 
as 

f{a) = A{a) + gpB{a) + glC{a) + glD{a) + g^Eia) + 0{gl). 
Each coefficient function is given to 10-loops as follows: 

A = 2a((T), 

B = 3%), 

C = -^q(3)(fT) - -C((7), 

D = fc(3)cic(3)(<7)-^ci,(5)(<7) + ^ci(<7), 

189C(3) 7C(3)^ , . 35C(5) 5C(7) 405 

^ = ec(3)(tx) + ^^e^(3)2((7) ec(5)(t^) + -^ec(7)(a)-— e(a), 

where 



r(i ^ _ ^ ^ ^ 7 _ 21 

715 \ 
V 



33 429 



65536(79 
1 / 1 

~ V ~ 2^ ' 8(j2 16(j3 ' 128(j4 256(j5 ' 1024(j6 2048a 



35 63 231 429 



_ 3/2^ 3 3 1 3 _ 3 _ 9 99 \ 

^^^^^ ~^ V^2^^8a2' 16^ ^ 128(7^ ~ 2560^ ^ 1024^ ~ 2048o^ ^ 32768ctV ' 
1/3 15 _ 35 315 _ 693 \ 

^ ~ ~ 2^^ 8^ "16^ ^128^" 256^7' 

, _ 1 1 1 5 7 21 \ 

_ 5/2/ 5 15 5 5 3 5 5 45 \ 

^^^^ ~ V ^ 2^ ^ 8^ ^ 16^ ~ 128^ ^ 256^ ~ 1024^6 ^ 2048(7^ ~ 32768(7V ' 
1/5 35 _ 105 \ 

2^^8^~T6^;' 
_ 1 / 1 3 5 35 \ 

""^^^^ " 7^1 2^ ^8^ "16^^128^7' 

5/2/ ^ 21 5 5 3 5_ 

ec(3)2 - (7 i + 2^ + 8^ + 16^ ~ 128^ ^ 256^ ~ 1024^ ^ 2048(7^ 
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2a 8(t2 



16^3 
35 



128(t4 256^5 
35 7 



1024(t6 
7 




+ 



16(t3 



+ 



128(t4 256cr5 



+ 



1024ct6 



2048ct7 



5 



+ 



32768CT' 



35 



) 



The coefficients, a, 6, c.^, and some of begin to appear at the lower loops (1,2 and 3 
loops). On the other hand, c, (i, e^(3) and e appear from 4, 6 and 8 loops and therefore we 
need higher order calculation to accumulate enough perturbative information. 
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Table Caption 

Table 1 The coefficient functions are evaluated via the stationarity condition. The 
number in () represents the loop level at which the evaluation is made. 

Figure Captions 

Fig.l Feynman diagrams dominant at each order in the small mass limit. 
Fig.2 Mdominantix) (3^ /{8C{3)) at 2-loop level is plotted for = 1 and g = 0.05. 



coefficient 










exact 


a 


1.12838(2) 


1.03027(4) 


1.01048(6) 


1.00417(8) 


1.00179(10) 1 


b 


0.75253(4) 


0.892268(6) 


0.948243(8) 


0.974006(10) 


1 


CC(3) 


0.56419(3) 


0.96975(5) 


0.99304(7) 


0.997867(9) 


1 


c 


0.300901(6) 


0.549637(8) 


0.718568(10) 




1 


C^C(3) 


1.12838(4) 


1.03027(6) 


1.01048(8) 


1.00417(10) 


1 




1.06455(5) 


1.00847(7) 


1.00192(9) 




1 


d 


0.08597(8) 


0.25135(10) 






1 


ec{3) 


0.75253(7) 


0.892268(9) 






1 


ec(3)2 


1.06455(6) 


1.00847(8) 


1.00192(10) 




1 


ec(5) 


0.56419(5) 


0.96975(7) 


0.99304(9) 




1 


ec(7) 


0.423142(3) 


0.619066(5) 


0.983706(7) 


0.997431(9) 


1 


e 


0.0191(10) 








1 



Table 1 



16 



